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Up to now, transverse quantum effects (usually labelled as "quantum imaging" efïects) which 
are generated by nonlinear devices inserted in resonant optical cavities have been calculated using 
the "thin crystal approximation" , i.e. taking into account the effect of diffraction only inside the 
empty part of the cavity, and neglecting its effect in the nonlinear propagation inside the nonlinear 
crystal. We introduce in the present paper a theoretical method which is not restricted by this 
approximation. It allows us in particular to treat configurations closer to the actual experimental 
ones, where the crystal length is comparable to the Rayleigh length of the cavity mode. We use 
this method in the case of the confocal OPO, where the thin crystal approximation predicts perfect 
squeezing on any area of the transverse plane, whatever its size and shape. We find that there 
exists in this case a "coherence length" which gives the minimum size of a detector on which perfect 
squeezing can be observed, and which gives therefore a limit to the improvement of optical resolution 
that can be obtained using such devices. 

PACS numbers: 42.50.Dv, 42.65.Yj, 42.60.Da 



I. INTRODUCTION 

Nonlinear optical elements inserted in optical cavities 
have been known for a long time to produce a great va- 
riety of interesting physical effects, taking advantage of 
the field enhancement effect and of the feedback provided 
by a resonant cavity 0,0- In particular, a great deal of 
attention has been devoted to cavity-assistcd nonlinear 
transverse effects, such as pattern formation Q and spa- 
tial soliton generation [5J- More recently the quantum 
aspccts of these phenomena have begun to be studied, 
mainly at the theoretical level, under the general name 
of " quantum imaging" , especially in planar or confocal 
cavities. 

Almost all the investigations relative to intra-cavity 
nonlinear effects, both at the classical and quantum level, 
have been performed within the mean field approxima- 
tion, in which one considers that the different interacting 
fields undergo only weak changes through their propaga- 
tion inside the cavity, in terms of their longitudinal and 
transverse parameters. This almost universal approach 
simplifies a great deal the theoretical investigations, and 
numerical simulations are generally needed if one wants 
to go beyond this approximation 9]. It implies in par- 
ticular that diffraction is assumed to be negligible inside 
the nonlinear medium, which limits the applicability of 
the method to nonlinear media whose length l c is much 
smaller than the Rayleigh length zr of the cavity modes 
z c (so called "thin" medium). This is a configuration 
that experimentalists do not like much: they prefer to 
operate in the case l c ~ z c which yields a much more effi- 
cient non linear interaction for a given pump power |lOj | . 



If one wants to predict results of experiments in realistic 
situations, one therefore needs to extend the theory be- 
yond the usual thin nonlinear medium approximation, 
and take into account diffraction effects oceurring to- 
gether with the nonlinear interaction inside the medium. 

The effects of simultaneous diffraction and nonlinear 
propagation have already been taken into account in 
the case of free propagation, i.e. without optical cavity 
around the nonlinear crystal, and they have been found 
to have a direct influence on the shape of the propagating 
beam Q. These effects have also been studied in detail 
at the quantum level in the parametric amplifier case , 
and recently for the soliton case 01 • I n contrast, they do 
not play a significant role when the nonlinear medium is 
inserted in an optical cavity with non degenerate trans- 
verse modes, which imposes the shape of the mode. But 
they are of paramount importance in the case of cavi- 
ties having degenerate transverse modes, such as a plane 
or confocal cavity, which do not impose the transverse 
structure of the interacting fields, and which are used to 
generate multimode quantum effects. 

Within the thin crystal approximation, i.e. taking into 
account diffraction effects only outside the crystal, strik- 
ing quantum properties have been predicted to oceur in 
a d egen erate OPO below threshold using a confocal cav- 
ity [3 El : this device generates quadrature squeezed 
light which is multimode in the transverse domain. It 
was shown in the case of a plane pump that the level of 
squeezing measured at the output of such an OPO nei- 
ther depends on the spatial profile of the local oscillator 
used to probe it, nor on the size of the detection region. 
This implies that a significant quantum noise reduction, 
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in principle tending to perfection when one approaches 
the oscillation thrcshold from below, can be observed in 
arbitrarily small portions of the down-converted beam. 
Therefore in this model there is no limitation in the trans- 
verse size of the domains in which the quantum noise 
is reduced when the OPO works in the exact confocal 
configuration. Such a multimode squeezed light appears 
thereby as a very promising tool to increase the resolu- 
tion in optical images beyond the wavelength limit. 

It is therefore very important to make a more realistic 
theoretical model of this system, which is no longer lim- 
ited by the thin crystal approximation, to see whcther 
the predicted local squeezing is still present in actual ex- 
perimental realizations in which the crystal length is of 
the order of the Rayleigh range of the resonator. This is 
the purpose of the present paper, in which we will show 
that the presence of a long crystal inside the resonator 
imposes a lower limit to the size of the regions in which 
squeezing can be measured ("coherence area"), which is 
proportional to WçI c /zr, where w c is the cavity beam 
waist, l c is the crystal length and zr the Rayleigh range 
of the resonator. 

The following section (section II) is devoted to the gen- 
eral description of the model that is used to treat the 
effect of diffraction inside the crystal, using the assump- 
tion that the single pass nonlinear interaction is weak in 
the crystal. We then describe in section III the method 
that is used to determine the squeezing spectra measured 
in well-defined homodyne detection schemes. We give in 
sections IV and V the results for such quantities respec- 
tively in the near field and in the far field, and conclude 
in section VI. 



II. THE MODEL 

A. Assumptions of the model 

Let us consider a confocal cavity, that for simplicity we 
take as a ring cavity of the kind shown schematically in 
Fig. n( 14] 15]). It is formed by four plane mirrors and 
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FIG. 1: Confocal ring cavity. The mirrors transmit the pump 
wave and reflect the signal wave, with the exception of mirror 
M c that partially transmits the signal. 



two lenses having a focal length equal to one quarter of 
the total cavity length, and symmetrically placed along 
the cavity, so that the focal points coincides at two posi- 
tions C and C. It contains a type I parametric medium 
of length l c , centered on the point C (see figure). It is 
pumped by a field A p of frequency 2uj s having a Gaussian 
shape and focused in the plane containing the point C. 
In such a plane the variation of the mean envelope with 
the transverse coordinate x is given by : 



A p {x) = A p exp(- J -^-) 



(1) 



We assume that the mirrors are totally transparent for 
the pump wave, and perfectly reflecting for the field at 
frequency u> s , except for the coupling mirror M c , which 
has a small transmission t at this frequency. The system 
was described in [Ï3| under the thin parametric medium 
approximation. We will follow here the same approach, 
generalized to the case of a thick parametric medium of 
length l c . The intracavity signal field at frequency lo s 
is described by a field envelope operator B(x,z), where 
z is the longitudinal coordinate along the cavity (z = 
corresponding to plane C), obeying the Standard equal 
time commutation relation at a given transverse plane at 
position z: 



[B(x, z, t), B T (x', z, t)] = S(x - x'). 



(2) 



As we are only interested in the regime below threshold 
and without pump dcplction, the pump field fluctuations 
do not play any role. 

In a confocal resonator the cavity resonances corre- 
sponds to complete sets of Gauss-Laguerre modes with 
a given parity for transverse coordinate inversion; we as- 
sume that a set of cavity even modes is tuned to reso- 
nance with the signal field, and that the odd modes are 
far off-resonance. It is then useful to introduce the even 
part of the field operator: 

B + (x, z, t) = \{B(x, z, t) + B(-x, z, i)] , (3) 
which obeys a modified commutation relation: 
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[B+(x, z, í), B\(x>, z, t)} = -[5(x - x') + 8(x + x')] (4) 

and can be written as an expansion over the even Gauss- 
Laguerre modes: 

È+(x,z,t) = ^2 f P ,i( x , z )àp,i(z,t), (5) 

p,leven 

where à Pt i(z, t) is the annihilation operator of a photon 
in mode (p, l) at the cavity position z and at time t. 

The interaction Hamiltonian of the system in the in- 
teraction picture is given by 

Hint = ^f I' dz 1 [ [ d 2 x'{A P (x\z')[Bl(x\z',t)] 2 -h.c.}, 



IJ2 



(6) 
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where g is the coupling constant proportional to the sec- 
ond order nonlinear susceptibility x ■ This equation 
generalizes the thin medium parametric Hamiltonian of 
Ref.Efi. 



Evolution equation in the image plane 
(near-fleld) 



B. 



In previous approaches 0, 0, the crystal was as- 
sumed to be thin, so that one could neglect the lon- 
gitudinal dependence of Ap and B + along the crystal 
length in the Hamiltonian ©. This cannot be done in 
a thick crystal. We will nevertheless make a simplify- 
ing assumption which turns out to be very realistic in 
the c.w. regime, with pump powers below 1W. We as- 
sume that the nonlinear interaction is very weak, so that 
it does not affect much the field amplitudes in a single 
pass through the crystal. We will therefore remove the 
z dependence of the operators à Pt i in Eqs.|0 |SJ, assum- 
ing à Pt i(z, t) = à p ,i(z = 0, t) = à p j(t), where z = is the 
crystal/cavity center C. The longitudinal variation of the 
signal operator B is then only due to diffraction and is 
described by the well-known z dependence of the modal 
functions / p ,;(x, z). This assumption leads to a rather 
simple expression of the commutator for the B + field at 
different positions inside the crystal: 

[È+(x,z,t),Èl(x\z',t')] = G;(z-z';x,x'). (7) 

Here G+(z;x, x') is the symmetrized part of the Fresnel 
propagator G(z;x, x'), describing the field linear propa- 
gation inside the crystal: 



G+ (z; x, x') = - [G(z; x, x') + G(z; x, -x')] . (8) 



with 



jb , X X' fl x| 



(9) 



where k s = n s uj s /c is the field wavenumber, with n s being 
the index of refraction at frequeney u) s , and we have intro- 
duced a walk-off term, present only if the signal wave is 
an extraordinary one, described by the two-dimensional 
walk-off angle p s . 

It is now possible to derive the time evolution of the 
field operator B(x, z, t) due to the parametric interaction. 
We will for example calculate it at the mid-point plane 
z = of the crystal : 



dB+ 



(x,0,t) 



= 9 



J J dVif mí (x,x").B^(x",0,í) 

(10) 



with the integral kernel Ki nt given by : 

x iní (x,x») = y J C/ f dz ' J J dV 



(11) 



In the limit of a thin crystal considered in Refs.[TÜ IÏ3|. 
Eq. I|ll·|l is replaced by the simpler expression : 



(x,0,í)| mí = gAp(x)Bl(x,0,t) 



(12) 



In the thin crystal case (Eq. I|12|l ). the parametric inter- 
action is local, i.e. the operators at different positions 
of the transverse plane are not coupled to each other, 
whereas in the thick crystal case ÍEa. (|l·lf> ). the paramet- 
ric interaction mixes the operators at different points of 
the transverse plane, over areas of finite extension. Note 
however that operators corresponding to different z vàl- 
ues are not coupled to each other, because of our as- 
sumption of weak parametric interaction. This situation 
is very close to the one considered in refs. 00,13] for para- 
metric down-conversion and amplification in a single-pass 
crystal, where finite transverse coherence areas for the 
spatial quantum cffects arise because of the finite spatial 
emission bandwidth of the crystal. In a similar way, in 
our case the spatial extension of the kernel Ki nt will turn 
out to give the minimum size in which spatial correlation 
or local squeezing can be observed in such a system. The 
analogy will become more evident in the next section, 
where we will explicitly solve the propagation equation 
of the Fourier spatial modes along the crystal. 

In order to get the complete evolution equation for the 
signal beam, one must add the free Hamiltonian evolution 
of the intracavity beam and the damping effects. This 
part of the treatment is Standard |17|. and is identical 
to the case of a thin crystal inserted in a confocal cavity 
[Ï3| . The final evolution equation reads: 



dB^ 
~dt 



(x,0,í) = - 7 (l + iA)B + (x,0,í) 



(13) 



dVüf int (x, x»)flt( x ",0, í) + ^/2 7 B +i „(x, 0, t) 



where 7 is the cavity escape rate, A the normalized cavity 
detuning of the even family of modes closest to resonance 
with the signal field, and Bi n the input field operator. 

In order to evaluate the coupling kernel, let us first take 
into account the diffraction of the pump field, focussed 
at the center of the crystal, z = 0. It is described by the 
Fresnel propagator G p (z;x, x'), equal to (JÜJ) when one 
replaces k s by the pump wavenumber k p , and the signal 
walk-off angle p s with the pump walk-off angle p p . One 
then gets : 



if íní (x,x") = 



rle/2 




/ dz' 


[ dV J 







d 2 y 



(14) 



A p (y)G p (z\ 0; x\ y)G* + (z', 0; x', x)G* + (z', 0; x\ x") 

Assuming for simplicity exact collinear phase matching 
k p = 2fc s , and neglecting the walk-off of the extraordi- 
nary wave , four of the five integrations can be exactly 
performed, and one finally gets: 



A P (x', z')G*Az'; x', x)G* («'; x', x") . 



if íní (x,x") 



A p (^±^)A(x-x») 



4 



+^( ZL V^) A ( X + X ") 



with 



A(x±x") 



dz> , 



!f |X±X'T 



(15) 



(16) 



-/c/2 



It can be easily shown that the function A(x±x") tends 
to the usual two-dimensional distribution <5(x±x") when 
l c — > 0, and that it can be written in terms of the integral 
sine function Si(x) = /* [3 



A(x±x") 



2tïL 



2 - H 2l c ) 



(17) 



This expression shows us that A takes negligible vàlues 
when |x± x"| > J Fig[21plots A as a function of 
the distance Ix ± x" I scaled to 



C. Evolution equation in the spatial Fourier 
domain (far-fleld) 

In this section we will investigate the intracavity dy- 
namics of the spatial Fourier amplitude of the signal field, 
which will offer an alternative formulation of the prob- 
lem. Fourier modes can be observed in the far-field plane 
with respect to the crystal center C, which in turn can 
be dctcctcd in the focal plane of a Iens placed outside the 
cavity. Let us introduce the spatial Fourier transform of 
the signal field envelope operator 

B + (q,M) = I ^fl + (x,z,t)e-«l· x 

= 1 \È(q, Z ,t) + È(-q,z,t)] (20) 



Equation (|14|) becomes: 



7171. 



wc 



lc 



n s z c 



(18) 



where wc and zq are the cavity waist and Rayleigh range, 
respectively. This expression shows that when the crys- 
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FIG. 2: Evaluation of the coupling kernel. A given by Eq. 
1161 1 is plotted as a function of |x±x" | scaled to the coherence 
length CHI. The first zero of A is obtained for the value 1.37 
of the coordinate. 

tal length is on the order of the Rayleigh range of the 
resonator, the transverse coherence length is on the or- 
der of the cavity waist. Recalling that the pump field 
has a Gaussian shape of waist w p , in order to have a 
multimode operation one must therefore use a defocussed 
pump, with w p ^> w c , or alternatively use a crystal much 
shorter than the Rayleigh range of the resonator, which is 
detrimental for the oscillation threshold of the OPO. The 
relevant sealing parameter of our problem is therefore 



l coh 



2n s 



lc 



(19) 



where z p is the Rayleigh or diffraction length of the pump 
beam. This parameter sets the number of spatial modes 
that can be independently excited, and it will turn out to 
give also the number of modes that can be independently 
squeezed. 



~dt 



■(q,0,i) = -7(l + *A)B+(q J 0,í) + 



(21) 



g J d^"K mí (q,q")^(q",0,í) + v/2 7 S +m (q,0,í), 

where the coupling Kernel Ki nt (q, q" ) is the Fourier 
transform of the kernel \ib\ with respect to both argu- 
ments. Straightforward calculations show that 



/v„„(q.q') = ïï (A p (q + q')sinc[^-|S_^|2] 



+À p (q - q')sinc[^_ 1 5L+S! j 2]^ (22) 
where A p is the spatial Fourier transform of the Gaussian 

~ 2 2 

pump profile QJ, i.e. À p (q) = ^fA p exp (-|q| 2 ^). 

The result (1221) can be also derived by solving the prop- 
agation equation of the pump and signal wave inside a 



X 



(2) 



crystal directly in the Fourier domain and in the 
limit of weak parametric gain. We will follow here the 
same approach as in and an d write the propa- 
gation equation in terms of the spatio-temporal Fourier 
transform field operators Aj (q, cj, z) of the pump (J = p) 
and signal (J — s) waves. Since the cavity linewidth is 
smaller by several orders of magnitude than the typical 
frequeney bandwidth of the crystal, the cavity filters a 
very small frequeney bandwidth around the carrier fre- 
queney u s of the signal; moreover, we have assumed that 
the pump is monochromatic, so that we can safely ne- 
glect the frequeney argument in the propagation equa- 
tions, which take the form 

^ (q, z) = ik jz (q)A, (q, z) + íf L (q, z) , (23) 

where is the nonlinear term, arising from the second 
order nonlinear susceptibility of the crystal. kj Z (q) — 

q 2 is the projection along the z-axis of the 
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wavevector, with kj = kj(u>j,Q) being the wave-number, 
which for extraordinary waves depends also on the prop- 
agation direction (identified by q). For the pump wave, 
we assume an intense coherent beam, that we suppose 
undepleted by the parametric down-conversion process 
in a single pass through the crystal, so that 



i p (q, z) A p (q, z) = e* k ^ z A p (q,0) 



(24) 



where we take the crystal center as the reference plane 
2 = 0. For the signal, the propagation equation is more 
easily solved by setting A s (q, z) — exp (ik sz (q)z)à s (q, z). 
The evolution along z of the operator à s is only due to the 
parametric interaction and is governed by the equation 
(see e.g.0 and 0] for more details): 

^|(q,*) = f / dV^(q + q',0)at(q',z) e W)^, 

(25) 

where u/l c is the parametric gain per unit length, and 
we have introduced the phase mismatch function 



5(q,q') = fc pz (q + q') - fc sz (q) - k sz (q'), 
Equation H25|) has the formal solution 



(26) 



as(q '"2 )+ L 



dz' I d 2 q' 



A p (q + q',0)àt(q',z')e^qq') z '. (27) 



Assuming a weak parametric efficiency er -C 1, we can 
solve this equation iteratively. At first order in a the 
solution reads: 

a,( q ,|) = a,(q,-|) + «7 J A'tfi(q,q')4(q',o), 

(28) 



with 



^i(q, q') = À p (q + q', 0)sinc 



(29) 



We observe that in the paraxial approximation kj Z (q) w 

2 

kj — pj ■ q — , where pj is the walk-off angle and kj = 
rijUJj/c. The phase mismatch function is hence given by: 



%,q') 



k p - 2k s + (p s - p p ) ■ (q + q') 



|q + q' 



(30) 



Assuming exact phase matching k p — 2k s , and neglecting 
the walk-off term, the argument of the sine function in 
Eq. (|2*ül) becomes 



%,q')| 



(31) 



In this way we start to recover the result of the Hamilto- 
nian formalism used to derive Eas. l|22|l and l|14|) . where, 



however, the effect of walk-off and phase mismatch were 
neglected for simplicity. Indeed, it is not difficult to 
show that the variation of the intracavity field operator 
-B+(q, 0, t) per cavity round trip time r, due to the para- 
metric interaction in a single pass through the crystal 



: AB+(q,0,í) 



-ür 1 (q ) -q , )]SÍ.(q',0,í)<32) 



This approach permits us to understand the physical ori- 
gin of the sine terms in the coupling kernel of Eq. I|22|) 
(which are the Fourier transform of the A terms in Eq. 
H151 ). that is the limited phase-matching bandwidth of 
the nonlinear crystal. For a crystal of negligible length, 
phase matching is irrelevant and there is no limitation in 
the spatial bandwidth of down-converted modes, whereas 
for a finite crystal the cone of parametric fluorescence 
has an aperture limited to a bandwidth of transverse 
wavevectors Aq « í/l C oh ^ l/V-Mc- As a consequence 
of the confocal geometry, the cavity ideally transmits all 
the Fourier modes, so that the only limitation in spatial 
bandwidth is that arising from phase matching along the 
crystal. 

We notice that if the pump is defocussed enough, the 
phase-matching limitation results in a limitation of the 
spot size oc 1/lcoh in the far- field with respect to the cav- 
ity center. Inside this spot, modes are coupled because 
of the finite size of the pump beam (the terms oc Ap in 
Eq. (1221) ) , inside a region of size oc w~ 1 . The relevant pa- 
rameter which sets the number of Fourier modes that can 
be independently excited is again given by b = w^H 2 coh 
(see Eq. (JTÜJ)). 



III. HOMODYNE DETECTION AND 
SQUEEZING SPECTRUM 

A. Homodyne detection scheme in the far field and 
near field 

The method used for measuring the noise-spectrum 
outside the cavity is a balanced homodyne detection 
scheme [2jj. We will use two configurations: the near- 
field configuration (x-position basis described in II. B) 
and the far-field configuration (q-vector basis described 
in II. C). The complete detection scheme in the near-field 
case is schematically shown in Fig. |3J The two matching 
lenses of focal length / image the crystal/cavity center 
plane C onto the detection planes D and D'. The image 
focal plane F of the first Iens coincides with the object 
focal plane of the second one, and represents the far-field 
plane with respect to the cavity center C. In planes C,F,D 
the signal field has its minimum waist, and it has a flat 
wave front. 

The detection scheme in the far field is obtained by 
using only one Iens as depicted in Fig. The focal 



6 



D|" O 

f „ f £ f „ f "> "" 




FIG. 3: Balanced homodyne detection scheme in the near 
field. Two matching lenses of focal f are used to image the 
cavity center C at the detection planes D and D' 

lcngth / Iens is used to image the far field plane with 
respect of the cavity center C onto the detection plane 
D. 




FIG. 4: Balanced homodyne detection scheme in the far field. 
A matching Iens of focal f is used to make the far field image 
of the cavity center C at the detection planes D and D' 

The symmetrical beam-splitter BS (reflection and 
transmission coefhcients r = -k= and t = -4=) mixes the 

V2 V2' 

output signal field with an intcnsc stationary and coher- 
ent beam ol(x, z), called local oscillator (LO). Note that 
all the fields being evaluated at the beam-splitter loca- 
tion, we will omit the z-dependence in the following. The 
difference photocurrent is a measure of the quadrature 
operator: 

En(0)= í d X [B out (x : n)a* L (x)+B out+ (x,-iï)a L ( X )l33) 

J det 

where det is the reciprocal image of the photodetection 



region at the beamsplitter plane, and assumed to be iden- 
tical for the two photodetectors. We have also assumed 
here that the quantum efficiency of the photodetector is 
equal to 1. Here B out is the sum of its odd and even part: 

B out (x, fi) = B° + at (x, fi) + B mt (x, fi) (34) 

The fluctuations 5Eh(ÇI) of the homodyne field around 
steady state are characterized by a noise spectrum: 

/+oo 
dn'{SE H (Q)SE H (íl')) = N + S(íl) (35) 
-oo 

where Eh is normalized so that N gives the mean photon 
number measured by the detector 

N= i dx\a L {x)\ 2 (36) 

J det 

N represents the shot-noise level, and S is the normally 
ordered part of the fluctuation spectrum, which accounts 
for the excess or decrease of noise with respect to the 
Standard quantum level. 



B. Input/output relation 

The relation linking the outgoing fields B^ lt (yi,t) 
with the intrac avity and input fields at the cavity in- 
put/output port[Ï3 is: 

B%* (x, í) = ^B ± (x, í) - B^ (x, t) (37) 

Equation (13) in the near-field (or (21) in the far field 
case) is easily solved in the frequency domain, by intro- 
ducing: 

5™ /ouí (x,fi) = / -^B" /mt (x,í)e- ,!!i 

Taking into account the boundary condition (37), we ob- 
tain the input/output relation: 



[ífi + 7 (l + ïA)][B7 í (x,fi)+ J B™(x,fi)] = 2 7 5™(x,fi) + 
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d 2 x K int (x,x ) 2 7j B™+(x , -fi) 



ifi + 7(1 - iA) 
d 2 x" jKlt (x , x" ) (Bf (x" , fi) + 57* (x" , fi) ) 

I 



(38) 



In the case of a thin crystal in the near field[l3j or a 
plane pump in the far field, this relation describes an 
infinite set of independent optical parametric oscillators. 
In these cases the squeezing spectrum can be calculated 



analytically as we will see in the following. But in other 
cases, this relation links all points in the transverse plane. 
In order to get the input/output relation, we have to 
inverse relation (38) by using a numerical method. 
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C. Numerical method 

In order to inverse relation (38) by numerical means, 
we need to discretise the transverse plane in order to 
replace integrals by discrete sums. For the sake of sim- 
plicity, we will only describe here the solution in the sin- 
gle transverse dimension model: the cavity is assumed 
to consist of cylindrical mirrors, so that the the trans- 
verse fields depend on a single parameter, y. In this 
case the electromagnètic fields are represented by vectors 
and the interaction terms by matrices. Straightforward 
calculations show that we can introduce the interaction 
functions U(y,y') and V(y,y') (calculated at resonance 
A = and at zero frequency in near-field or far-field con- 
figurations ) linking two different points in the transverse 
plane, so that relation (38) becomes: 



B° + Ut (y) 



dy'U(y,y')Bf(y') 



dy'V(y,y')B™+(y>) (39) 



Since we assumed that the odd part of the output field 
is in the vacuum state, B™ 1 gives no contribution to the 
normally ordered part of the spectrum S, which can be 
calculated by using the input /output relation (21) for 
the even part of the field, and by using the commutation 
rules for the even part: 



[Bf /out (x,t),B^ /out+ (x',t')} = 
hs(x-x')±S(x + x')]S(t-t') 



(40) 



In the following, we will assume, as in Refs. 01 > that 
the local oscillator has a constant phase profile ^l(x) = 
ifL, so that oil(x) = \a L (x)\e l· ' PL . We obtain the ordered 
part of the spectrum, normalized to the shot noise: 



S(0) = 1 
N J det dy\a L \ 



dxdx' 



det 2 



dy\a L (x)\\a L (x')\ [{V(x, y)V(x', y) + V(x, y)V(x, -y)) 

+ cos(2<p L )(U(x, y)V(x', y) + U(x, y)V(x', -y)) (41) 



Now, knowing the U(y,y') and V(y,y') interaction func- 
tions, we are able to calculate the squeezing spectrum in 
both near and far-field cases. 



IV. SQUEEZING SPECTRUM IN THE 
NEAR-FIELD 

In this section, we use the near-field homodyne detec- 
tion (Fig. |2l described in Ref.0]. As already said in 
part II, in the near field, the thick crystal couples pixels 
contained in a region whose size is in the order of l co h 
(18). 

Let us study first the case of a plane-wave pump and 
a plane-wave local oscillator. As pointed out in ref[Ï2|. 
in this case and in the thin crystal approximation, the 
level of squeezing does not depend on the width of the 
detection region. Fig. [S] shows results predicted for a 
measurement performed with a circular detector of radius 
Ap centered on the cavity axis (which is a symmetric 
detection pointed out in [ï^)- We represent the 

squeezing spectrum at zero frequency as a function of 
the size of the detector, sealed to the coherence length 

lcoh = \j ^r- We can see that for Ap < l c , the squeezing 
tends to zero when Ap — ► 0, as already predicted. For 
larger vàlues of the detector size, perfect squeezing can 
be achieved. We can also see that the squeezing evolution 
is comparable to the A function evolution (Fig. |2J). 



In the more realistic case of finite size pump, the 

w 2 

squeezing level depends on the parameter b = p- 2 - = 

coh 

2n s | £ , as pointed out in part I. Fig. [ü] represents the 
squeezing spectrum at zero-frequeney as a function the 
detector radius, normalized to l co h, for different b param- 
eters, using a plane local oscillator. As already seen in 
Fig. for Ap — > 0, the noise reduction effect tends 
to zero. But we see now that there is also no squeezing 
effect for large vàlues of the detector radius, because of 
the finite size of the pump, as already shown in |13|. 

Fig. [7|shows theoretical results in the case of a detector 
consisting of two symmetric pixels (pixel of size equal to 
the coherence length), for different b vàlues, in function 
of the distance between the two pixels. For large vàlues 
of p, the noise level goes back to shot noise because of the 
finite size of the pump, as already depicted in reference 
[13|. But now, for small p vàlues, the squeezing does not 
tend to zero, as in the thin crystal case. 



V. SQUEEZING SPECTRUM IN THE 
FAR-FIELD 

In this section, we will consider the spatial squeezing 
spectrum in the far field (Fig. and in the q-vector 
basis. As already said in section II. C, the coupling be- 
tween q- vectors modes is now due to the finite length of 
the pump. We will see that a new coherence length l co hf 



8 



0.6 




coli 



FIG. 5: Squeezing spectrum at zero-frequency, normalized to 
the shot noise, as a function of the detector radius (scaled to 

Icah) • 
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FIG. 7: Squeezing spectrum at zero-frequency, normalised to 
the shot noise, as a function of the pixel distance between the 
two pixels p from the cavity axis (scaled to l C oh), plotted for 
several vàlues of b 



0.8 

Z 0.6 
— 

o 

ll 

G 

> 0.4 
0.2 




Ap/1 



coh 



FIG. 6: Squeezing spectrum at zero-frequency, normalised to 
the shot noise, as a function of the radial amplitude of the 
detector scaled to l coh , plotted for several vàlues of b 



appears in such a case, given by: 



Icohf oc — 



. We will successively investigate two configurations: the 
plane wave pump regime (where the squeezing spectrum 
can be calculated analytically) , and the case of a finite 
pump size (where a numerical method is necessary). 



A. Plane wave pump regime in the far fleld 

In order to evaluate the far field case, we introduce 
the spatial Fourier transforms of the electromagnètic field 
temporal frequency components: 



d 2 x 



i/out i 



In the case of a plane wave pump, A p (x, z) = A p , so that 
equation(22) becomes: 



dt 



■(q,0,t) = - 7 [(l + iA)B + (q,0,í) + v/2 7j B +m (q,0,í) 



A pS mc(|^)^(q,0,í)] 



(42) 



This equation, which does not mix different q vàlues, can 
be solved analytically. It is similar to equation (14) in 
reference |l3j . Taking into account the boundary condi- 
tion 

B° ± ut (q, í) = y^B ± (q, í) - Éf (q, t) (43) 
We obtain: 

Sf*(q,0) = tf(q,Q)S;«(q,Q) + V(q, íï)B*?+(-q, -íï) 

(44) 

where 

[1 - »(A - íï/ 7 )][l - »(A + 0/ 7 )] + ^ S znc 2 (^ff ) 



E/(q,í2) 



[1 + i(A + n/ 7 )][l - i(A - í2/ 7 )] - Ajsinc 2 ^) 

(45) 
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and 



V(q,Ü) 



(46) 



2A p sinc( l -j^) 



[1 + i(A + fi/ 7 )][l - i(A - fi/ 7 )] - A2 S mc 2 (^f ) 

In the case of plane wave regime, the input/output re- 
lation in the spatial Fourier space describes therefore an 
infinitc set of independent optical parametric oscillators 
below threshold. This can be simply understood: the q- 
vector basis is the eigenbasis of the diffraction, so that 
no coupling between q-vector modes due to the crystal 
appears. 

Let us now consider the homodyne-detection scheme, 
schematically shown in (Fig. 0J. The Iens provides a 
spatial Fourier transform of the output field B ou t(x, fi), 
so that at the location of plane D the field B® ut (x, fi) is: 



2tt - 



B» ut ( x ,n) = —B out (— x ,n) 



,2tt 



(47) 



In this plane, B® ut (x, fi) is mixed with an intense station- 
ary and coherent beam a£ (x) = jj&Loi^xf , wnere 
«i(x) has a gaussian shape, with a waist w^o- The ho- 
modyne field has thus an expression similar to the near 
field case, where functions of x are now replaced by their 
spatial Fourier transforms: 



E H (fi) 



(48) 



dq[B° ut (q, fi)à* LO (q) + B° ut+ (q, -íl)à LO (q)} 



del 



This analogy shows that, in the case of a local oscillator 
that has an even parity with respect to coordinate inver- 
sion, the squeezing spectrum is given by (like in [Ï3^: 



smc(^— ) < are quite outside the phase matching 
curve, so that the choice fLo{q) = § everywhere should 
give good results. The squeezing spectrum at resonanec 
and zero frequeney, for í^loÍq) = § can be analytically 
calculated and, as a function of the radius r of a detector 
centered on the optical axis is given by: 



V(r,0) 
N 



where 



rr/ro / - 

Jo wexp(- 

rr/ro 

* / uexp 
Jo 



ro 



(52) 



k,u z 




1 + A p sinc(u 2 ) 
1 — A p sinc(u 2 ) 



(53) 



Fig. |H1 shows the results obtained in the case of two dif- 
ferent detection configurations: the V curve shows re- 
sults in the case of a circular detector of variable radius r 
(sealed to ro), using a local oscillator waist wlo = tq. As 
already said in part //, the limitation of the squeezing 
level is due to the non perfect phase matching along the 
crystal. For r > r , the squeezing level decreases. So, in 
the plane wave pump regime in the far field, the thickness 
of the crystal has a role comparable with the finite size 
of the pump in the near field, as reported in [l3T |. The 
R curve shows results obtained in the case of two small 
symmetrical pixels and a plane wave local oscillator as 
a function of the pixel distance from the cavity axis r, 
sealed to ro. We can see that the noise level goes back to 
the shot noise level for r > ro, because to the non perfect 
phase matching along the crystal. 



V(Ü) 



dq{\à L o(q)\ 2 {l-a(q)}} 



de / 



dq{\à LO (q)\ 2 <r(q)R(q,m (49) 



del 



where the noise spatial density R(q, fi) is given by: 

R{q,n) = |t/(q,fi) + e 2í ^°( í «y*(q,-fi)| 2 (50) 
and where 



<7(q)= / dq'í+(q,q') 

Jdet 



(51) 



In order to minimize R(q, fi), the local oscillator phase 
should be chosen as fLo(q) = ar9 ^ u ^'^ v ^' n ^ . In par- 
ticular, at resonance and at zero frequeney Í7(q, 0) and 
V(q, 0) are real and the optimal local oscillator phase 

would correspond to (fLo(q) = f , when smc(-^S-) > 

0, and (fiLo(q) = 0, when sinc(^l·) < 0, which is 
not indeed very practical. However, modes for which 



B. Squeezing spectrum in the far field case and 
finite size pump regime 

When one takes into account the finite size of the 
pump, a coupling between different q vectors appear, and 
one needs to solve equations numerically, as in the near 
field case. A new coherence length l co hf appears in the 
far field: l cohf = 

Fig. |ÏÏ|shows the evolution of the squeezing spectrum at 
zero frequeney, and at resonance, for different b param- 
eters, in function of the detector radius sealed to lcohf 
. We see the same evolution as in the analytical case, 
except that the noise level tends to shot noise for small 
vàlues of the detector. 

Fig. ^3 shows the results obtained in the case of two 
symmetrical pixels (pixel of size equal to the coherence 
length lcohf), for different b vàlues, in function of the dis- 
tance between the two pixels p. The evolution is similar 
to the one given by fig(6) for large distances. But there 
is here also a decrease of the squeezing effect for small 
distances. 
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FIG. 8: Squeezing spectrum normalized to the shot noise,at 
zero-frequency, at resonance, in the plane pump regime and 
far field case, for two measurement configurations. V is ob- 
tained using a circular detector of radial amplitude r (scaled 
to ro). R is obtained using a pair of symmetrical pixels in 
function of the pixel distance from the axis r (scaled to ro) 
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FIG. 10: Squeezing spectrum normalized to the shot noise 
at zero-frequency, and at resonance, as a function of the dis- 
tance between the two pixels p (scaled to the coherence area 
Icohf), in the finite pump regime and far field approach and 
for different vàlues of b. 
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FIG. 9: Squeezing spectrum normalized to the shot noise at 
zero-frequency, and at resonance, as a function of the radial 
amplitude of the detector Ap(scaled to the coherence area 
Icohf), in the finite pump regime and far field approach and 
for different vàlues of b. 



VI. DISCUSSIONS AND CONCLUSIONS 

We have seen that when one takes into account the 
effect of diffraction inside the nonlinear crystal in a con- 
focal OPO, the local squeezing predicted for any shape 
and size of the detectors in the thin crystal approximation 
is now restricted to detection areas lying within a given 
range, characterized by a coherence length l co h- This 



prediction introduces serious limitations to the success 
of an experiment, and must be taken into account when 
designing the experimental set-up. With the purpose of 
producing a light beam that is squeezed in several ele- 
mentary portions of its transverse cross-section, either a 
crystal short comparcd to zr should be used or, alter- 
natively, a defocussed pump, with a waist much larger 
than the cavity waist. In both cases the efficiency of the 
non linear coupling is reduced. For instance, with lcm 
long crystal, l co h is equal to 40/ito, and one must choose 
a pump waist much larger than this value in order to ob- 
serve multimode squeezing (the number of modes being 

w 2 

roughly equal to the ratio b = jt 2 -). This defocussed 

l coh 

pump will imply a much higher threshold for the OPO 
oscillation, which is multiplied by a factor also close to b. 
The conclusion of this analysis is that one cannot have 
multimode squeezing "for free", and that with a given 
pump power, one will be able to to excite a number of 
modes which is roughly equal to the ration of the injected 
pump power to the threshold power for single mode op- 
eration. 
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